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Abstract 

Spectrum sensing, i.e., detecting the presence of primary 
users in a licensed spectrum, is a fundamental problem 
in cognitive radio. Since the statistical covariances of re- 
ceived signal and noise are usually different, they can be 
used to differentiate the case where the primary user's 
signal is present from the case where there is only noise. 
In this paper, spectrum sensing algorithms are proposed 
based on the sample covariance matrix calculated from 
a limited number of received signal samples. Two test 
statistics are then extracted from the sample covariance 
matrix. A decision on the signal presence is made by com- 
paring the two test statistics. Theoretical analysis for the 
proposed algorithms is given. Detection probability and 
associated threshold are found based on statistical theory. 
The methods do not need any information of the signal, 
the channel and noise power a priori. Also, no synchro- 
nization is needed. Simulations based on narrowband sig- 
nals, captured digital television (DTV) signals and multi- 
ple antenna signals are presented to verify the methods. 



1 Introduction 

Conventional fixed spectrum allocation policy leads to low 
spectrum usage in many of the frequency bands. Cogni- 
tive radio, first proposed in (1J, is a promising technology 
to exploit the under-utilized spectrum in an opportunistic 
manner [31 [3J 01 [5] . One application of cognitive radio is 
spectral reuse, which allows secondary networks/users to 
use the spectrum allocated/licensed to the primary users 
when they are not active [6]. To do so, the secondary 
users are required to frequently perform spectrum sens- 
ing, i.e., detecting the presence of the primary users. If 
the primary users are detected to be inactive, the sec- 
ondary users can use the spectrum for communications. 
On the other hand, whenever the primary users become 
active, the secondary users have to detect the presence 
of those users in high probability, and vacate the chan- 
nel within certain amount of time. One communication 
system using the spectrum reuse concept is IEEE 802.22 
wireless regional area networks (WRAN) [7], which oper- 
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ates on the VHF/UHF bands that are currently allocated 
for TV broadcasting services and other services such as 
wireless microphone. Cognitive radio is also an emerging 
technology for vehicular devices. For example, in [5J cog- 
nitive radio is proposed for underwater vehicles to fully 
use the limited underwater acoustic bandwidth, and in [9] 
it is used for autonomous vehicular communications. 

Spectrum sensing is a fundamental task for cognitive 
radio. However, there are several factors which make 
spectrum sensing practically challenging. First, the sig- 
nal to noise ratio (SNR) of the primary users may be 
very low. For example, the wireless microphones oper- 
ating in TV bands only transmit signals with a power of 
about 50mW and a bandwidth of 200 kHz. If the sec- 
ondary users are several hundred meters away from the 
microphone devices, the received SNR may be well below 
— 20dB. Secondly, multipath fading and time dispersion 
of the wireless channels make the sensing problem more 
difficult. Multipath fading may cause the signal power 
fluctuates as large as 20 — 30dB. On the other hand, co- 
herent detection may not be possible when the time dis- 
persed channel is unknown, especially when the primary 
users are legacy systems which do not cooperate with the 
secondary users. Thirdly, the noise/interference level may 
change with time, which yields noise uncertainty. There 
are two types of noise uncertainty: receiver device noise 
uncertainty and environment noise uncertainty. The re- 
ceiver device noise uncertainty comes from [TQl [TTJ [12] : (a) 
non-linearity of components; and (b) time- varying ther- 
mal noise in the components. The environment noise un- 
certainty may be caused by transmissions of other users, 
either unintentionally or intentionally. Because of noise 
uncertainty, in practice, it is very difficult to obtain the 
accurate noise power. 

There have been several sensing methods, including the 
likelihood ratio test (LRT) [13] , energy detection method 
[HI M [El El H2l [15], matched filtering (MF)-based 
method [TT], [121 HS1 US] and cyclostationary detection 
method [17l [TBI E] > each of which has different require- 
ments and advantages/disadvantages. Although LRT is 
proved to be optimal, it is very difficult to use it in prac- 
tice, because it requires exact channel information, and 
distributions of source signal and noise. In order to use 
LRT for detection, we need to obtain the channels, and sig- 
nal and noise distributions first, which are practically in- 
tractable. MF-based method requires perfect knowledge of 
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the channel responses from the primary user to the receiver 
and accurate synchronization (otherwise its performance 
will be reduced dramatically) [151 116] . As mentioned ear- 
lier, this may not be possible if the primary users do not 
cooperate with the secondary users. Cyclostationary de- 
tection method needs to know the cyclic frequencies of the 
primary users, which may not be realistic for many of the 
spectrum reuse applications. Furthermore, this method 
demands excessive analog to digital converter (ADC) re- 
quirement and signal processing capabilities Energy 
detection, unlike the other two methods, does not need 
any information of the signal to be detected and is robust 
to unknown dispersed channel and fading. However, en- 
ergy detection requires perfect knowledge of noise power. 
Wrong estimation of the noise power leads to SNR wall 
and high probability of false alarm pH EJ [H Q21 [20] . As 
pointed out earlier, the estimated noise power could be 
quite inaccurate due to noise uncertainty. Thus, the main 
drawback for the energy detection is its sensitiveness to 
noise uncertainty [TQl [TT1 [12l [15] . Furthermore, while en- 
ergy detection is optimal for detecting independent and 
identically distributed (iid) signal [13], it is not optimal 
for detecting correlated signal, which is the case for most 
practical applications. 

In this paper, to overcome the shortcoming of energy de- 
tection, we propose new methods based on statistical co- 
variances or auto-correlations of the received signal. The 
statistical covariance matrices or auto-correlations of sig- 
nal and noise are generally different. Thus this difference 
is used in the proposed methods to differentiate the signal 
component from background noise. In practice, there are 
only limited number of signal samples. Hence, the detec- 
tion methods are based on the sample covariance matrix. 
The steps of the proposed methods are as follows. First, 
the sample covariance matrix of the received signal is com- 
puted based on received signal samples. Then two test 
statistics are extracted from the sample covariance matrix. 
Finally, a decision on the presence of the signal is made by 
comparing the ratio of two test statistics with a threshold. 
Theoretical analysis for the proposed algorithms is given. 
Detection probability and associated threshold for deci- 
sion are found based on statistical theory. The methods 
do not need any information of the signal, the channel and 
noise power a priori. Also, no synchronization is needed. 
Simulations based on narrowband signals, captured digi- 
tal television (DTV) signals and multiple antenna signals 
are presented to evaluate the performance of the proposed 
methods. 

The rest of the paper is organized as follows. The detec- 
tion algorithms and theoretical analysis are presented in 
Section II. Section III gives the performance analysis and 
finds thresholds for the algorithms. Theoretical compar- 
ison with the energy detection is also discussed in this 
section. Simulation results for various types of signals 
are given in Section IV. Conclusions are drawn in Sec- 
tion V. Finally, a pre-whitening technique is given in the 
Appendix. 



Some notations are used as follows: boldface letters 
are used to denote matrices and vectors, superscript (-) T 
stands for transpose, \ q denotes the identity matrix of or- 
der q, E[-] stands for expectation operation. 

2 Covariance Based Detections 

Let x c (t) = s c (t) + rj c {t) be the continuous-time received 
signal, where s c (t) is the possible primary user's signal 
and T] c (t) is the noise. rj c {t) is assumed to be a station- 
ary process satisfying E(?y c (i)) = 0, E(?7^(i)) = and 
F,(r) c (t)i] c (t + t)) = for any r ^ 0. Assume that we are 
interested in the frequency band with central frequency 
f c and bandwidth W. We sample the received signal at 
a sampling rate f s , where f s > W. Let T s = l// s be 
the sampling period. For notation simplicity, we define 
x(n) = x c (nT s ), s(n) = s c (nT s ) and r/(n) = rj c (nT s ). 
There are two hypothesizes: Tlo, the signal does not exist; 
and Tii, the signal exists. The received signal samples un- 
der the two hypothesizes are given respectively as follows 

P33 Kg H21 ng : 

Ho : x(n) = r](n), (1) 
Hi: x(n)=s(n)+T)(n), (2) 

where s(n) is the transmitted signal samples passed 
through a wireless channel consisting of path loss, mul- 
tipath fading and time dispersion effects, and r](n) is 
the white noise which is independent and identically dis- 
tributed (iid), and with mean zero and variance az. Note 
that s(n) can be the superposition of the received signals 
from multiple primary users. No synchronization is needed 
here. 

Two probabilities are of interest for spectrum sensing: 
probability of detection, Pd, which defines, at the hypoth- 
esis Tii, the probability of the sensing algorithm having 
detected the presence of the primary signal; and proba- 
bility of false alarm, Pf a , which defines, at the hypothesis 
Tio, the probability of the sensing algorithm claiming the 
presence of the primary signal. 

2.1 Covariance Absolute Value (CAV) 
Detection 

Let us consider L consecutive samples and define the fol- 
lowing vectors: 

x(n) = [ x(n) x(n-l) ■■■ x(n - L + 1) ] T (3) 

s(n) = [ s(n) s(n-l) ••• s(n - L + 1) ] T ,(4) 

77(71) = [ r){n) T](n — I) ■■■ rj(n — L + 1) ] T .(5) 

The parameter L is called smoothing factor in the follow- 
ing. Considering the statistical covariance matrices of the 
signal and noise defined as 

K x = E[x(n)x T (n)], (6) 
R 5 = E[s(n)s T (n)], (7) 
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we can verify that 



(8) 



If the signal s(n) is not present, R s = 0. Hence the off- 
diagonal elements of R x are all zeros. If there is signal 
and the signal samples are correlated, R s is not a diagonal 
matrix. Hence, some of the off-diagonal elements of R^ 
should be non-zeros. Denote r nm as the element of matrix 
R-r at the nth row and mth column, and let 
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(9) 
(10) 



Then, if there is no signal, T\jT 2 = 1. If the signal is 
present, Ti/T 2 > 1. Hence, the ratio T 1( /T 2 can be used 
to detect the presence of the signal. 

In practice, the statistical covariance matrix can only 
be calculated using a limited number of signal samples. 
Define the sample auto-correlations of the received signal 
as 
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A W = ^T E x(m)x(m-l), J = 0,1,---,L-1, (11) 



m— 



where iV s is the number of available samples. The statis- 
tical covariance matrix R^ can be approximated by the 
sample covariance matrix defined as 



A(0) 
A(l) 



A(l) 
A(0) 



X(L-l) A(L-2) 



A(L-l) 
X(L - 2) 

A(0) 



.(12) 



Note that the sample covariance matrix is symmetric and 
Toeplitz. Based on the sample covariance matrix, we pro- 
pose the following signal detection method. 

Algorithm 1 The covariance absolute value (CAV) de- 
tection algorithm 

Step 1. Sample the received signal as described above. 

Step 2. Choose a smoothing factor L and a threshold 
7i , where 71 should be chosen to meet the requirement for 
the probability of false alarm. This will be discussed in the 
next section. 

Step 3. Compute the auto- correlations of the received 



signal X(l), I = 0, 1, 
variance matrix. 
Step 4- Compute 

Ti(JV a ) 
T 2 (N S ) 



, L — 1, and form the sample co- 



L L 

r nm {N s )\, 

n—l m— 1 
L 

1 nn 



L L 

= EE 



1 L 

= iE 



(13) 
(14) 



where r nm (N s ) are the elements of the sample covariance 
matrix TL X (N S ). 

Step 5. Determine the presence of the signal based 
on Ti(N s ), T 2 (N S ) and the threshold 71, i.e., if 
Ti(N s ) /T 2 (N S ) > 71, signal exists; otherwise, signal 
does not exist. 



Remark. The statistics in the algorithm can be cal- 
culated directly from the auto-correlations A(Z). How- 
ever, for better understanding and easing the mathemati- 
cal derivation for the pre-whitcning later in Appendix A, 
here we choose to use the covariance matrix expression. 

2.2 Theoretical Analysis for CAV Algo- 
rithm 

The proposed method only uses the received signal sam- 
ples. It does not need any information of the signal, the 
channel and noise power as a priori. Also, no synchroniza- 
tion is needed. 

The validity of the proposed CAV algorithm relies on 
the assumption that the signal samples are correlated, that 
is, R s is not a diagonal matrix (some of the off-diagonal 
elements of R s should be non-zeros). Obviously, if the 
signal samples s(n) are iid, then R s = (t^Il- In this case, 
the assumption is invalid and the algorithm cannot detect 
the presence of the signal. 

However, usually the signal samples should be corre- 
lated due to the following reasons. 

(1) The signal is oversampled. Let T be the Nyquist 
sampling period of the signal s c (t) and s c (nT ) be the 
sampled signal based on the Nyquist sampling rate. Based 
on the sampling theorem, the signal s c (t) can be expressed 
as 



3 (t)= J2 s c (nT a )g(t-nT a ), 



(15) 



where g(t) is an interpolation function. Hence, the sig- 
nal samples s(n) = s c (nT s ) are only related to s c (nT ). 
If the sampling rate at the receiver f s > 1/Tq, that is, 
T s < To, then s(n) = s c {nT s ) must be correlated. An 
example of this is the narrowband signal such as the wire- 
less microphone signal. In a 6 MHz bandwidth TV band, a 
wireless microphone signal only occupies about 200 KHz. 
When we sample the received signal with sampling rate 
not lower than 6 MHz, the wireless microphone signal is 
actually over-sampled and therefore highly correlated. 

(2) The propagation channel has time dispersion, thus 
the actually signal component at the receiver is given by 



=(*) 



j — t 



h{r)so(t - r)dr, 



(16) 



71=1 



where s (t) is the original transmitted signal and h(t) is 
the response of the time dispersive channel. Since the 
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sampling period T s is usually very small, the integration 
(|16p can be approximated as 



s c (t) w T s ^2 h(kT s )s (t - kT a 



(17) 



Hence, 



s c (nT s ) w T s 



E 

fc=-R"o 



h(kT s )s ((n - k)T s 



(18) 



where [ifo^aj ^iT s ] is the support of the channel response 
h(t), that is, /i(t) = for t g [K T s ,KiT s ]. For time dis- 
persive channel, K\ > K Q , thus the signal samples s(nT s ) 
are correlated even if the original signal samples So(nT s ) 
could be iid. 

(3) The original signal is correlated. In this case, even if 
the channel is a flat fading channel and no oversampling, 
the received signal samples are correlated. 

Another assumption for the algorithm is that the noise 
samples are iid. This is usually true if no filtering is used. 
However, if a narrowband filter is used at the receiver, 
sometimes the noise samples will be correlated. To deal 
with this case, we need to pre-whiten the noise samples or 
pre-transform the covariance matrix. A method is given 
in Appendix A to solve this problem. 

The computational complexity of the algorithm is as 
follows. Computing the auto-correlations of the received 
signal requires about LN S multiplications and additions. 
Computing T\{N a ) and T 2 (-/V s ) requires about L 2 addi- 
tions. Therefore, the total number of multiplications and 
additions is about LN S + L 2 . 

2.3 Generalized Covariance Based Algo- 
rithms 

Based on the same principle as CAV, generalized covari- 
ance based methods can be designed to detect the signal. 
Let ip\ and ip2 be two non-negative functions with multiple 
variables. Assume that 

V>i(a) > 0, for a f 0; ^(0) = 0; 

^a(b) >0, forb^O; ^ 2 (0) = 0. 

Then, the following method can be used for signal detec- 
tion. 

Algorithm 2 Generalized Covariance based Detection 

Step 1. Sample the received signal as described above. 

Step 2. Choose a smoothing factor L and a threshold 
72, where 72 should be chosen to meet the requirement for 
the probability of false alarm. 

Step 3. Compute the sample covariance matrix R X (N S ). 

Step 4- Compute 



T±(N S 
T 3 (N S 



tp2{r nn (N s ), n 
T 4 (A S )+Vi (r„ 



= 1,- 



n =/= m) 



(19) 
(20) 



Step 5. Determine the presence of the signal based 
on T^(N S ), T4(N S ) and the threshold 72. That is, if 
T3(N s )/Ti(N s ) > 72, signal exists; otherwise, signal does 
not exist. 

Obviously the CAV algorithm is a special case of the 
generalized method when ipi and ip2 are absolute sum- 
mation functions. As another example, we can choose 
ipi (a) = a T a and ^(b) = b T b. For this choice, 
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n—1 in— 1 

1 L 
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(21) 
(22) 



2.4 



Spectrum sensing Using Multiple An- 
tennas 



Multiple antenna systems have been widely used to in- 
crease the channel capacity or improve the transmission 
reliability in wireless communications. In the following, 
we assume that there are M > 1 antennas at the receiver, 
and exploit the received signals from these antennas for 
spectrum sensing. In this case, the received signal at an- 
tenna i is given by 

H.q : Xi(n) = r)i(ri), (23) 
Hi: Xi{n) = Sl {n) + r ll {n). (24) 

In hypothesis Hi, Sj(n) is the signal component received 

by antenna i. Since all Si(n) are generated from the same 

source signal, thus Si(n) are correlated for i. It is assumed 

that ?7i(n)'s are iid for n and i. 

Let us combine all the signals from the M antennas and 
define the following vectors: 

x(n) = [ xi (n) ■■■ XM(n) xi(n — 1) ■■• xu(n — 1) 

■•■ xi(n-L + l)--- x M (n-L + l) ] T , (25) 
s(n) = [ si(n) ■•■ s M (n) si (n—1) ■■■ s M (n-l) 

■•• si{n- L + l)--- s M (n-L + l) } T , (26) 

77(71) = [ r?i(n) ■•■ r) M (n) 7/1 (n — 1) ■■• r?j|f(n— 1) 

... rj^n- L + l)--- VM (n-L + l) ] T . (27) 

Note that equations ^ to ([5]) are a special case (M = 1) 
of the above equations. Defining the statistical covariance 
matrices using the same way as those in (J6)) and ([7]), we 
obtain 



(28) 



Except the different matrix dimensions, the equation 
above is the same as equation (jHJ. Hence, the CAV algo- 
rithm and generalized covariance based method described 
above can be directly used for multiple antenna case. 

Let So(n) be the source signal. The received signal at 
antenna i is 



JV, 

E 

fe=0 



= E h i{k)s {n ~k) + Vi(n), i = 1,2, • • • ,M, (29) 
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where hi(k) is the channel responses from the source user 3.1 Statistics Computation 

to antenna i at the receiver. Define 



h(n) = [hx(n), h 2 (n), ■ 
h(0) h(N) ■ 







h(0) 



, h M (n)} 


h(N) 



(30) 
(31) 



where N = max(A^) and hi(n) is zero-padded if Ni < N. 

i 

Note that the dimension of M is ML x (TV + L). We have 



— MR, 



(32) 



where R So = E(§o§^) is the statistical covariance matrix 
of the source signal, where 



so 



[ s (n) s (n - 1) 



s (n-N-L + l) ] 



(33) 

Note that the received signals at different antennas are 
correlated. Hence, using multiple antennas increase the 
correlations among the signal samples at the receiver and 
make the algorithms valid at all cases. In fact, at the 
worst case when all the channels are flat-fading, that is, 
Ni = N2 = ■ ■ ■ = Nm = 0, and the source signal sample 
so(n) is iid, we have R s = ct^HH t , where H is a ML x L 
matrix defined above. Obviously, R s is not a diagonal 
matrix and the algorithms can work. 



Performance Analysis 
Threshold Determination 



and 



For a good detection algorithm, Pd should be high and 
Pf a should be low. The choice of the threshold 7 is a 
compromise between the Pd and Pf a . Since we have no 
information on the signal (actually we even do not know 
if there is signal or not), it is difficult to set the threshold 
based on the Pd. Hence, usually we choose the threshold 
based on the Pf a . The steps are as follows. First we set 
a value for Pf a . Then we find a threshold 7 to meet the 
required Pf a . To find the threshold based on the required 
Pfa, we can use either theoretical derivation or computer 
simulation. If simulation is used to find the threshold, we 
can generate white Gaussian noises as the input (no sig- 
nal) and adjust the threshold to meet the Pf a requirement. 
Note that the threshold here is related to the number of 
samples used for computing the sample auto-correlations 
and the smoothing factor L, but not related to the noise 
power. If theoretical derivation is used, we need to find the 
statistical distribution of Ti(N s )/T 2 (N s ), which is gener- 
ally a difficult task. In this section, using central limit 
theorem, we will find the approximations for the distribu- 
tion of this random variable and provide the theoretical 
estimations for the two probabilities, Pd, Pfa, as well as 
the threshold associated with these probabilities. 



Based on the symmetric property of the covariance matrix, 
we can rewrite Ti(N s ) and T 2 (N S ) in ([13j) and iJHJ) as 



L-l 



T X (7V S ) = A(0) + z ]>>-Z)|A(0|, (34) 
1=1 

T 2 {N S ) = A(0). (35) 



Define 



X ; = [ x(N s -l-l) ••• x(-l) Y , (36) 
Si = [s(N s -l-l) ■■■ s(-l) f, (37) 

m = [ n {N s -i-i) ■■■ v (-i) f. (38) 

Let the normalized correlation among the signal samples 
be 

ai =E[s(n)s(n- l)]/a 2 , (39) 

where a 2 is the signal power, a 2 = E[s 2 (rt)]. |a/| defines 
the correlation strength among the signal samples, here 
^ lad ^ !• Based on the notations, we have 



1 

W. 
1 



A(/! = —XTx l = ±(sT + r,Z)(S l •„,.! 



Obviously, 

E(A(0)) = a 2 + <j 2 , 
E(A(0) = a t a 2 s , I = 1, 2, • • • , L - 1. 
Now we need to find the variance of X(l). Since 



(40) 



(41) 
(42) 



A2 W = jp ( s o s < + s o m + nSSt + vlvi) 2 



1 



(SfS,) +(S^) +(^S ; ) + (r,^) 

+2 (S^S,) (S^)+2(S^S ; ) (ttfS,) 
+2 (S^Sz) (ttffi,)+2(Sfo) (t|2"S,) 
+2 (Sfo,) (r, T om ) + 2 (?7o S/) (itfij,)] , (43) 
it can be verified that 



E((S^) 2 ) 
E((^ ) 2 ) 
E((^) 2 ) 
E((SjS0(S^)) 



E((S r S )( J7o r f7o )) 
EftSjTSzXfjjfo)) 
E ((Sjfo)(,#S,)) 



E^S,) 2 ) =N s a 2 a 2 v , (44) 
(N 2 + 2N s )a% (45) 
N s a% l = l,--- ,L-1, (46) 
E((SjS ; )(rjJSO) 
E((S^)(r?^)) 
E (fao SOfafoi)) =0, (47) 
N 2 a 2 a 2 , (48) 
0, I = 1,2,- • ■ ,L- 1, (49) 
a 2l (N s ~ l)a 2 a 2 v (50) 
ma 2 . (51) 
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Based on these results, we can easily obtain the following Proof. Based on the central limit theorem, we have 
two lemmas. 



E(|A(/)|) 



1 



Lemma 1 When there is no signal, we have 



+ 00 



1 



2 . = — == / |vA;M + 9i|e 2 du 

E(A(0)) = a*, Var(A(0)) = — a% (52) J-c 



+ OC 



-oo 

-e,/V57 



E(A(i)) = 0, Var(A(/)) = — ^, 1 = 1, 1. (53) = J-o. ( — V — 0;)e - ~ 3 ~dti 
Lemma 2 When there is signal, we have + ,— / (v A;u + 0;)e 2 dit 



Var(A(0)) = Var (^S ) + ^ (2a 2 + o*) (55) 



2Q l f° 



E(A(i)) = a Jf r s 2 , (56) V2ir J-e,/^ 



e 2 du 



Var(A(0) = Varf-U^S,) = A /^£ ( 2 _ e -&\ + ^gi /° e -^# dw 



2(N S - l)a 2l 

-1- 

« = !,•■• ,L- 1 



7T 



2tt J-ej^Ki 



— K + 2<x s + a s (57) = 2 _ 



7T 



+ OO 



+|0,| - e~-du .(62) 

Note that Var f^S^S,] , / = 0,1,- •• depends on \ V tt J^i/v^T j 

signal properties. For large N s and low SNR, 

For simplicity, we denote E(A(Z)) by &i and Var(A(Z)) by 222 

A/. Note that usually N s is very large. Based on central ^ ^ (°"s + (7 ?)) |QjJ ^ 

limit theorem, A(^) can be approximated by the Gaussian N. 

distribution. c , , . 

So, we obtain (|60p 



Lemma 3 W-Tierc i/ie signa/ is no£ present, we have wg ]3 tam equation 



When there is no signal, 9; = and A; = j^-at- Hence, 



Theorem 1 When there is no signal, we have 



E (l A (0l) = \/^< 1 = 1,2,-.. ,L-1. (58) 



When the signal is present, we have 



E(T x (iV s )) = (i + (L-l)yZjL)a* (63) 



E(|A(0l) = ^{2-e-^ E(T 2 (7V S )) = o*, (64) 

+m (l-Jir --*d«)(B9) = ^ 

When there is signal, and for large N s , we have 

E(Ti(N s )) w erf 



77 7|e t |/VSr 
/ = 1,2,--. ,Z-1. 



for Zarge iV s and Zow SNR, 

E(|A(/)|) « W^-^ + a^h-e^ " '-' V 



2~~ / -n 2 \ 1 L 

•"■jVs V / 0^2 1 ^ i-i 



|6z| I 1 - \/ - / e-Vdu , (60) u 1=1 



7T 

Z = 1,2,... ,L- 1 



(66) 

I E(T 2 (N s ))=a 2 s +at (67) 



where 

|a / |SNR % /iV7 OAT „ o- 2 



Var(T 2 (7V s) ) = V ar (^S ) + ^ (2a 2 + a 2 ) . (68) 



77 = 1 " — — — -, SNR = — . (61) For large N s , Ti(N s ) and T 2 (N S ) approach to Gaussian 

1 + SNR a-2 distributions. 
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Proof. Equations (|63|) to 168)) are direct results from 
Lemma [T] to [3] Noting that is a summation of N s 
random variables, when N s is large, based on the central 
limit theorem it can be approximated by Gaussian distri- 
butions. From the definition of Ti(N s ) and T2(N S ), we 
know that they also approach to Gaussian distributions. 



□ 



3.2 Detection Probability and the Associ- 
ated Threshold 

From the theorem above, we have 



9rr 2 

lim E(T 1 (Ay)=a s 2 + ( x2 + -^( J L-0h|. (69) 



1=1 



For simplicity, we denote 



L-l 



at 



(70) 



i=i 



which is the overall correlation strength among the con- 
secutive L samples. When there is no signal, we have 



T!(iv s )/r 2 (iv s ) 



E(T!(iV s ))/E(r 2 (JV fl )) 



(71) 



Note that this ratio approaches to 1 as iV s approaches to 
infinite. Also note that the ratio is not related to the noise 
power (variance). On the other hand, when there is signal 
(signal plus noise case), we have 



T l {N s )/T 2 {N s ) a E(T!(Ar s ))/E(T 2 (iV s )) 



= 1 



SNR , 
SNR+ 1 



(72) 
(73) 



Here the ratio approaches to a number larger than 1 as 
N s approaches to infinite. The number is determined by 
the correlation strength among the signal samples and the 
SNR. Hence, for any fixed SNR, if there are sufficiently 
large number of samples, we can always differentiate if 
there is signal or not based on the ratio. 

However, in practice we have only limited number of 
samples. So, we need to evaluate the performance at fixed 
N s . First we analyze the Pf a at hypothesis Hq. The 



probability of false alarm for the CAV algorithm is 

P fa = P(T 1 (iV s )>7 1 r 2 (iV s )) 



= P[T 2 {N S )<—T 1 (N S ) 
V 7i 



P \T 2 {N S ) < — ( 1 + (L — 1) 



N s ir I a " 



= P 



- f 



_^_(j2 



where 



Q(*) 



1 



-|-oc 



~ u2 / 2 du. 



(74) 



V2tt Jt 

For a given Pf a , the associated threshold should be chosen 
such that 



^( 1 + ( i - 1 ) V / E) 



f 



That is, 



71 = 



-Q _1 ( P /a) 



(75) 



(76) 



Note that here the threshold is not related to noise 
power and SNR. After the threshold is set, we now cal- 
culate the probability of detection at various SNR. For 
the given threshold 71, when signal presents, 



P, 



P{T 1 {N s )> ll T 2 {N s )) 
P (t 2 (N s ) < ^-Ji(N s ) 



P T 2 (N S ) < — E(Ti(JV a )) 

V 71 



P 



T 2 (N S ) - <7g - a 2 v ^E(Ti(iV s )) - a 2 - a 2 ' 



< 



1-Q 



y/Vfa(T 2 (N s )) ^/Vm(T 2 {N s )) 
/iE(T x (iV s ))-(7?-^\ 



v/Var(T 2 (7V s )) 
For very large N s and low SNR, 
2cr 2 

Var(T 2 (Ag)«^(2a 2 +< 7 2 ) 



(77) 



2(<x s 2 +a 2 ) 2 



Hence, we have a further approximation 



l-Q 



1-Q 



71 

y/Wa 

J_ , Tj,SNR 
71 "f" 7i(SNR+l) 



(78) 
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Obviously, the Pd increases with the number of samples, 
N s , the SNR and the correlation strength among the signal 
samples. Note that 71 is also related to N s as shown above, 
and lim 71 = 1. Hence, for fixed SNR, Pd approaches 

N s ^oo 

to 1 when N s approaches to infinite. 

For a target pair of Pd and Pf a , based on ([78]) and ([76]) . 
we can find the required number of samples as 



(T L SNR) 2 



For fixed Pd and Pf a , N c is only related to the smoothing 
factor L and the overall correlation strength T^. Hence, 
the best smoothing factor is 



L best = mm{7V c }, 



(80) 



which is related to the correlation strength among the sig- 
nal samples. 

3.3 Comparison with Energy Detection 

Energy detection is the basic sensing method, which 
was first proposed in 14] and further studied by others 
[T01 HU 021 [15] . It does not need any information of the 
signal to be detected and is robust to unknown dispersive 
channel. Energy detection compares the average power 
of the received signal with the noise power to make a de- 
cision. To guarantee a reliable detection, the threshold 
must be set according to the noise power and the number 
of samples [TTJ [10] [12] ■ On the other hand, the proposed 
methods do not rely on the noise power to set the thresh- 
old (see equation ([75]) ). while keeping other advantages of 
the energy detection. 

Accurate knowledge on the noise power is then the key 
of the energy detection. Unfortunately, in practice, the 
noise uncertainty always presents. Due to the noise un- 
certainty []JJ El [H], the estimated (or assumed) noise 
power may be different from the real noise power. Let the 
estimated noise power be cr 2 , - 
uncertainty factor (in dB) as 



•ticf^. We define the noise 



B = max{10 log 10 a}. 



(81) 



It is assumed that a (in dB) is evenly distributed in an in- 
terval [— B, B] [11]. In practice, the noise uncertainty fac- 
tor of a receiving device is normally 1 to 2 dB [TTJ [20]. The 
environment/interference noise uncertainty can be much 
higher [llj . When there is noise uncertainty, the energy 
detection is not effective jTO] [TT] [T^] [2D]. Simulation in 
the next section also shows that the proposed method is 
much better than the energy detection when noise uncer- 
tainty is present. Hence, here we only compare the pro- 
posed method with the ideal energy detection (without 
noise uncertainty). 

To compare the performances of the two methods, first 
we need a criterion. By properly choosing the thresholds, 
both methods can achieve any given Pd and Pf a > if 



sufficiently large number of samples are available. The key 
point is how many samples are needed to achieve the given 
Pd and Pf a > 0. Hence, we choose this as the criterion 
to compare the two algorithms. For energy detection, the 
required number of samples is approximately [llj 



N,= 



2 (Q-\Pfg) - Q-\Pd)Y 

SNR 2 



(82) 



(79) Comparing ([79]) and (|82|). if we want N c < N e , we need 



T L > 1 + 



L — 1 



V^(Q- 1 (^/a)-Q- 1 (^))' 



(83) 



For example, if Pd = 0.9 and Pf a — 0.1, we need 
T l > 1 + j-gj- In conclusion, if the signal samples are 



highly correlated such that ([83]) holds, the CAV is bet- 
ter than the ideal energy detection; otherwise, the ideal 
energy detection is better. 

In terms of the computational complexity, the energy 
detection needs about N s multiplications and additions. 
Hence, the computational complexity of the proposed 
methods is about L times that of the energy detection. 

4 Simulations and Discussions 

In this section, we will give some simulation results for 
three situations: narrowband signals, captured DTV sig- 
nals |21] and multiple antenna received signals. 

First, we simulate the probabilities of false alarm (Pf a ) 
because Pf a is not related to signal (at Tio, there is no 
signal at all). We set the target Pf a — 0.1, and choose 
L = 10 and N s = 50000. We then obtain the thresholds 
based on the Pf a , L and N s . The threshold for energy de- 
tection is given in [TTJ. Table[JJgives the simulation results 
for various cases, where and in the following "EG-x dB" 
means the energy detection with x-dB noise uncertainty. 
The Pf a for the proposed method and energy detection 
without noise uncertainty meet the target, but the Pf a 
for the energy detection with noise uncertainty (even as 
low as 0.5 dB) far exceeds the limit. This means that the 
energy detection is very unreliable in practical situations 
with noise uncertainty. 

Secondly, we fix the thresholds based on Pf a and sim- 
ulate the probability of detection (Pd) for various cases. 
We consider two signal types as follows. 

(1) Narrowband signals. FM modulated wireless mi- 
crophone signal is used here (soft speaker) [22]. The cen- 
tral frequency is f c = 100 MHz. The sampling rate at the 
receiver is 6 MHz (the same as the TV bandwidth in USA) . 
Figure [TJ gives the simulation results (the correspond- 
ing Pfa is shown in Table [lj. Note that "CAV-theo" 
means the theoretical results given in Section III.B. Due 
to some approximations, the theoretical results does not 
exactly match the simulated results. The CAV is better 
than ideal energy detection (without noise uncertainty), 
which verifies our assertion in Section III. C. The reason 
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method 


EG- 2 dB 


EG-1.5 dB 


EG-1 dB 


EG-0.5 dB 


EG-OdB 


CAV 


Pfa 


0.497 


0.496 


0.490 


0.470 


0.102 


0.099 



Table 1: Probabilities of false alarm 




SNR (dB) 

Figure 1: Probability of detection for wireless microphone 
signal: N s = 50000 



0.9 - ; 




0.01 0.05 0.1 0.15 0.2 0.25 0.3 

Probability of false alarm 



is that, as we point out in Section II. B, the source signal is 
narrowband, and therefore their samples are highly corre- 
lated. As shown in the figure, if there is noise uncertainty, 
the Pd of the energy detection is much worse than that of 
the proposed method. Figure [2] gives the Receiver Operat- 
ing Characteristics (ROC) curve (Pd versus Pf a ) at fixed 
SNR = — 20dB. The performances of the methods at dif- 
ferent sample sizes (sensing times) are given in Figure [3J 
It is clear that CAV is always better than the ideal energy 
detection. 

To test the impact of the smoothing factor, we fix 
SNR=-20dB, P fa = 0.01, and N s = 50000, and vary the 
smoothing factor L from 4 to 14. Figure 2] shows the re- 
sults for the Pd- We see that the Pd is not very sensitive 
to the smoothing factor for L > 8. Noting that smaller L 
means lower complexity, in practice, we can choose a rel- 
atively small L. However, it is very difficult to choose the 
best L because it is related to signal property (unknown) . 
Note that energy detection is not affected by L. 

(2) Captured DTV signals. The real DTV signals 
(field measurements) are collected at Washington D.C., 
USA. The data rate of the vestigial sideband (VSB) DTV 
signal is 10.762M samples/sec. The recorded DTV signals 
were sampled at 21.524476M samples/sec and down con- 
verted to a low central IF frequency of 5.381119 MHz (one 
fourth the sampling rate). The analog-to-digital conver- 
sion of the RF signal used a 10-bit or a 12-bit A/D. Each 
sample was encoded into a 2-byte word (signed intl6 with 
a two's complement format). The multipath channel and 
SNR of the received signal are unknown. In order to use 
the signals for simulating the algorithms at very low SNR, 



Figure 2: Pd versus Pf a for wireless microphone signal: 
N s = 50000, SNR = -20dB 




Sample size 



Figure 3: Pd versus sample size (N s ) for wireless micro- 
phone signal: P fa = 0.01, SNR = -20dB 
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Smoothing factor 



-22 -20 -18 -16 -14 -12 -10 -8 -6 -4 -2 
SNR (dB) 



Figure 4: Pd versus smoothing factor for wireless micro- 
phone signal: P fa = 0.01, N s = 50000, SNR = -20dB 




SNR (dB) 



Figure 5: Probability of detection for DTV signal WAS- 
051/35/01: N s = 50000 

we need to add white noises to obtain various SNR levels 
[23]. 

Figure Ogives the simulation results based on the DTV 
signal file WAS-051/35/01 (receiving antenna is outside 
and 20.29 miles from the DTV station; antenna height is 
30 feet) [21]. The corresponding Pf a is shown in Table 
[TJ If noise variance is exactly known (B = 0), energy 
detection is better than the proposed method. However, 
as discussed in |T01 [TTJ |T2], noise uncertainty is always 
present. Even if noise uncertainty is only 0.5dB, the Pd of 
energy detection is much worse than that of the proposed 
method. 

In summary, all the simulations above show that the 
proposed method works well without using information of 
the signal, the channel and noise power. The energy detec- 
tion are not reliable (low probability of detection and high 
probability of false alarm) when there is noise uncertainty. 

Thirdly, we simulate the proposed algorithms with mul- 



Figure 6: Probability of detection using multiple antennas: 
Pfa =0.1, N s — 25000, one source signal 

tiple antennas/receivers. We consider a system with 4 
receiving antennas. Assume that the antennas are well 
separated (separation larger than half wavelength) such 
that their channels are independent. This assumption is 
only for simplicity. In fact, the proposed algorithms per- 
form better if the channels are correlated. Assume that 
each multipath channel hi(k) has 5 taps (iVj = 4) and all 
the channel taps are independent with equal power. The 
channel taps are generated as Gaussian random numbers 
and different for different Monte Carlo realizations. The 
source signal So(n) are iid and BPSK modulated. The re- 
ceived signal at antenna i is defined in (|29|) . The smooth- 
ing factor is L — 8 and the number of samples at each 
antenna is N s = 25000. We fix the Pf a — 0.1 at all cases. 
Figure [6] gives the Pd for three cases. From the figure we 
see that, when only one antenna's signal is used (M = 1), 
the method still works. This verifies our assertion in Sec- 
tion II. B that the method is valid even if the inputs are 
iid but the channel is dispersive. When 2 (M = 2) or 
4 (M = 4) antennas' signals are combined based on the 
method in Section II. D, the Pd is much better. The more 
the antennas we use, the better the Pd is. The optimal 
LRT [13] detection for M = 4 is also included as an upper 
bound for any detection methods. 

To check the performance of the methods for time vari- 
ant channels, we give a simulation result here. The time 
variant channel is generated based on the simplified Jake's 
model. Let fd be the normalized maximum Doppler fre- 
quency (DF). The time variant channel for simulation is 
defined as 

hi(n, k) = exp \ j2-Kn^-^f d \ h^k), k = 1, • • • , 5, (84) 

where hi{k) is the time invariant channel defined above. 
For different Doppler frequency from to 10 -2 , simulation 
result is shown in Figure [7] for M = 2. For fast time 
variant channels, the performance of proposed methods 
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-22 -20 -1 



SNR(dB) 



Figure 7: Probability of detection for time variant chan- 
nels: M = 2, P fa =0.1, N s = 25000 




SNR (dB) 



Figure 8: Probability of detection using multiple antennas: 
Pfa =0.1, N s = 25000, three source signals 



will degrade. 

We also simulate the situation when there are multiple 
source signals. Figure gives the results for the case of 
three source signals. Compared to Figure [6j here the re- 
sults do not change much. Hence, the proposed method is 
valid when there are multiple source signals. 

5 Conclusions 

In this paper, sensing algorithms based on the sample co- 
variance matrix of the received signal have been proposed. 
Statistical theories have been used to set the thresholds 
and obtain the probabilities of detection. The methods 
can be used for various signal detection applications with- 
out knowledge of the signal, the channel and noise power. 
Simulations based on the narrowband signals, captured 
DTV signals and multiple antenna signals have been car- 



ried out to evaluate the performance of the proposed meth- 
ods. It is shown that the proposed methods are in general 
better than the energy detector when noise uncertainty is 
present. Furthermore, when the received signals are highly 
correlated, the proposed method is better than the energy 
detector even the noise power is perfectly known. 

Appendix A 

At the receiving end, sometimes the received signal is fil- 
tered by a narrowband filter. Therefore, the noise embed- 
ded in the received signal is also filtered. Let r](n) be the 
noise samples before the filter, which are assumed to be 
iid. Let f(k), k = 0, 1, • • • , K, be the filter. After filtering, 
the noise samples turns to 



fj(n) 



K 

E 

k=0 



f(k)rj(n — k), n = 0, 1, ■ 



Consider L consecutive outputs and define 
fj( n ) = [fj(n), ■■■ ,fj(n-L + 1)] T . 



(85) 



(86) 



The statistical covariance matrix of the filtered noise be- 
comes 



R-T, 



E(fj(n)fj(n) T ) 



where F is a L x (L + K) matrix defined as 

~/(0) ••• f{K-l) /:7m ••• 







/(0) /(i) 



f(K) 



(87) 



(88) 



Let G = FF T . If analog filter or both analog and digital 
filters are used, the matrix G should be defined based on 
those filter properties. Note that G is a positive definite 
symmetric matrix. It can be decomposed to 



G = Q 2 , 



(89) 



where Q is also a positive definite symmetric matrix. 
Hence, we can transform the statistical covariance matrix 
into 



Q 



R, ; Q = a 2 v h 



(90) 



Note that Q is only related to the filter. This means that 
we can always transform the statistical covariance matrix 
R x in ^ (by using a matrix obtained from the filter) such 
that equation ([8]) holds when the noise has passed through 
a narrowband filter. Furthermore, since Q is not related 
to signal and noise, we can pre-compute its inverse Q 1 
and store it for later usage. 
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